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PHYSICAL PROPERTIFS OF STANDARD AGENTS, CANDIDATE AGENTS, AND 
RELATED COMPOUNDS AT SEVERAL TEMPERATURES (U) 


1. (U) INTRODUCTION 


(U) A problem facing users of physical chemical properties of chemical 
agents and related compounds has been the lack of a central source of reliable 
data. Data are often scattered in various technical reports, journals and 
laboratory notebooks. Many compilations of data fail to note original sources 
or important experimental details. As a result, the physical properties used 
in calculations in various technical reports are inconsistent and often not 
the best data available. 


(U) For over a decade, efforts of chemists in the Physical Organic 
Section have been directed toward providing an automated physical property 

data system for chemical agents and related compounds. The characteristics 

of the desired system are listed in table 1. To provide the maximum assistance 
to the user of physical chemical properties of chemical agents, every effort 
was made to obtain these characteristics. 


(U) This report briefly describes the physical property data system 
along with methods used to predict additional properties. A compilation of 
physical, chemical and thermodynamic data and correlations is provided for 
55 chemical agents and related compounds. A list of these compounds with 
their chemical structures is given in Appendix A. In Appendix B, the 
properties along with predicted properties are provided at temperatures of 
-40°, -209, 0°, 20°, 25° and 40°C. 
Table 1, (U) Characteristics of the Desir<d Data System 

Use automation to ease access and to allow for upgrading of input. 

Enter properties in formatted form so that each value is machine readable. 

-Evaluate all available data and use best data. 

Combine data to extend the experimental data range where appropriate. 

Calculate properties at any specified temperature. 


Provide easy to read printout with readily apparent units. 


Take all data from original sources with provided references. 


Provide data ranges and warnings where extrapolations are performed. 
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2. (Uv) PHYSICAL CHEMICAL DATA SYSTEM 
(u) All values for properties and equation parameters entered into the 
data system are formatted and identified by property code and compound number. 
Thus, all properties are machine readable except solubility which lists the 
solvent in the reference field. Table 2 lists the major components of an 80 
character input line. 
Table 2. (U) Components of a Data Input Line 
Compound code number 
Compound letter code 


Property or parameter code number 


Property or parameter value 


Temperature and/or upper or lower values 
of temperature range 


Reference field of up to 36 characters 


UNCLASSIFIED 


(u) Compounds are listed in order of their EA number. A system 
of numbered codes was established to identify individual properties and 
parameters, The compound numbers, names, and structures for compounds 
contained in this report are tabulated in Appendix A. 


(0) Considerable time and effort were spent locating, evaluating and 
selecting physical chemical properties to be included in this system. The 

data selected are considered the best available. When conflicting data were 
found, selection was made primarily on the basis of experimental procedure 

and reported sample purity. Therefore, the values in this system tend to be for 
purified rather than technical grade material. When properties had been 
extensively evaluated and documented in reviews, these values were used unless 
better data had become available. 


(U) Two types of equations are used in this system to describe physical 
properties versus temperature. Where possible, data sets covering different 
temperature ranges were combined to generate property parameters based on 
extended temperature ranges. 


(u) Antoine equations! are used to describe both vapor pressure and 
viscosity. The merits of the Antoine equation for fitting and extrapolating 
vapor pressure data are covered extensively in Appendix C. The Antoine equation 
has been found also to work well for fitting and interpolating viscosity data. 
Hovever, for many of the compounds, the Antoine constants for viscosity were 
generated from three points which covered a data range of only 25°C. Extensive 
extrapolation of viscosity for these compounds may result in a significant 
error. 
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(u) Density, refractive index, and surface tension are each expressed 
as a linear function of temperature. The method of least squares was used to 
generate equation parameters from the experimental data. 


(uv) The reference field allows up to 36 characters per property or set 
of parameters. Due to this space limitation, it was necessary to abbreviate 
references. These references should be available in the Chemical Systems 
Laboratory (CSL) Technical Library. Since some data are unpublished, it was 
necessary to list CSL notebook (NB) numbers. Several unreferenced data points 
are included in the system despite substantial efforts to reference all infor- 
mation. Additional information is sometimes provided in the reference field, 
e.g., sample purity and other characteristics. In cases where the property 
value was reported in the literature as a range, an average value was used and 
the range was reported in the reference field. 


(u) Numerous tests and warnings have been incorporated into this system 

for the following reasons: (1) to prevent the calculation of meaningless data, 
and (2) to warn the user when the calculated values fall outside the range of | 
the experimental data used in generating the input parameters, The most common 
warning is "THE ABOVE VALUES ARE EXTRAPOLATED OUT OF THE DATA TEMPERATURE RANGE." 
Due to the small temperature coefficients, extrapolations of the linear functions 
such as density and surface tension over broad temperature ranges are relatively 
reliable. The viscosity, and vapor pressure, along with properties derived from 
vapor pressure, undergo significant changes with temperature. Relatively poor 
values for these properties could result from large extrapolations, particularly 
when equations are derived from data collected over a narvrow temperature range. 
An extrapolation of vapor pressure over a 200° temperature range could result 

in errors of one magnitude or more in value. Tests and warnings also have 

been incorporated into the data system to prevent extrapolation of data from 

one physical state to another. Table 3 demonstrates how the system uses freezing 
points (fp) (or melting points) and boiling points (bp) to classify a specified 
temperature prior to calculating liquid properties. Appropriate statements 

and warnings are printed along with the calculated values. 


Table 3. (U) Test and Warnings 


Classification Directed Action 


bp < T Liquid properties are calculated at bp 
fp < T < bp Liquid properties are calculated at T 
(£p-259C) < T << fp Liquid properties are calculated at T and it is 


noted that these values are calculated for a 
supercooled liquid 


T < (£p-25°C) Liquid properties are calculated at the fp 
No fp Liquid properties calculated at T with warning 
that there is no fp available 
No bp Liquid properties calculated at T with warning 
that there is no bp available 24 
UNCLASSIFIED 
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3. (u) PREDICTIONS OF PHYSICAL PROPERTIES 


(u) Early in the development of the physical chemical property data 
system, it was realized that a great potential existed with this system for 
predicting additional properties. Bretsznajader's book, Prediction of Transport 
and Other Physical Properties of Fluids,? provided much help in selecting pro- 
cedures which required properties consistent with those contained within this 
data system. Various empirical relationships were added to the data system 

to allow the prediction of critical properties, the viscosity of vapor, and 

the diffusion coefficient of vapor into air. 


3.1 (U) Critical Properties. 


(U) The critical density, critical temperature, critical volume, and 
critical pressure are calculated using the method of Filippov.3 The property 
parameters used to calculate critical properties are the density equation 
constants A and B, and the molecular weight. The following equation allows the 
critical density (p.) in gm/ml to be calculated at T degrees Kelvin (°K). 


dp 
Pe * 0.253 \09 - T UF (1) 


where 9 is the density of compound at temperature T, and do/dT is the derivative 
of density with respect to T. 


Critical temperature (T,) in °K is calculated using the following 
relationship: 


= -1.95 de 
¢c do /aT X (2) 


Critical pressure (P,) in atmospheres (atm) is calculated as follows: 


RT cP ¢ 
3.83 M 


P. * (3) 


where R is 82.06 (ml atm/deg-mole), and M is the molecular weight. 


The critical volume (V,) is simply the gram molecular weight 
divided by the critical density. 


(u) Using Filippov's method, good estimates of critical properties should 
be obtained when the following criteria are met for the density equation parameters: 
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b. Experimental range is in an area where the vapor density is 
small compared to liquid density. 


4. 


-. 


c. dp/dT is constant over a broad temperature range. 


(u) For most liquids, do/dT is relatively constant. Highly associated 
liquids are an exception. Water exhibits an abnormality at 4°C where do/dT 

e changes signs. Above 49°C, the magnitude of do/dT increases significantly with 
increased temperature. 


(U) Filippov? found that errors in calculating critical properties 

usually did not exceed 2 percent for pop, and 3 percent for Tg. Literature 
values for critical properties of chemical agents contained in this report 

were not available; therefore, four test compounds were added to the system. 
Table 4 contains a listing of the calculated critical properties and literature 
values for these compounds. Input for the test compounds consisted of molecular 
weights and density equation parameters. Density parameters were derived from 
literature values of density at various temperatures over a range of 25 to 35°C. 
The 2-propanol was selected as a nonideal case. Due to hydrogen bonding, its 
density is a nonlinear function of temperature. 


(U) The estimated critical properties are in relatively good agreement with 
the literature values for the first three compounds. The estimated critical 
temperature for diethylamine is 22.5°C lower than the reported value; however, the 
difference in absolute temperature is only 5 percent. The estimated critical 
temperature for 2-propanol which is about 15 percent lower than the experimenral 
value can be attributed primarily to the nonlinearity of its density versus 
temperature relationship. 
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3.2 (U) Diffusion Coefficients. 


The method of Othmer and Chen’ is used to calculate the diffusion 
coefficient of a vapor into air. 


( 
4 
‘ 
: 
\ 
{ 
‘ 


Dig? = (2.52 K 107) w2*7* (4) 
a 
where Di» = the diffusion coefficient of component one into component 
? 2 
two in cm’ /sec 
M, = molecular weight of component one (vapor) 
M, = molecular weight of component two (air) 
P = pressure in atm 
Ve, = critical volume in ml of component one (vapor) 
Veo = critical volume in ml of component two (air) 
yp = viscosity of air in centipoise (cP) at the desired temperature. 
(U) The V, of component one is calculated as described in the pre- 
ceding section on critical properties. The value used for the molecular 
weight of air is 28.958 while the Pe Of air is reported to be 0.35 g/m1.9 
Thus, the V, of air is calculated to be 82.71 ml. In this report, all cal- 
culations of diffusion coefficients are made at one atm. 
(U) The viscosity of air is calculated at the required temperatures 
using the Sutherland10 equation for the influence of temperature on the 
viscosity of gases. When the constants for air* are substituted into the 
equation, the following relationship is obtained: 
Mage 0.001488 12/2 (5) 
7 T + 120 | 
where Voir is the viscosity of air at T in cP. 
9 (U) Using diffusion coefficient literature values for 50 systems at 


25°C, Othmer and Chen’ found an average absolute deviation of 5.75 percent. 
In addition, they compared their calculated values for an air-water system 
between 25 and 1220°C to observed values and found an average error of 7.61 
percent. Citing experimental difficulties, they concluded that their method, 
as well as several more complex methods, should give values within experi- 
mental error. 
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(uU) Diffusion coefficients for the vapor into air for the four test 
compounds were estimated in the physical property data system and included 

in table 4. Estimated values agree within 7 percent for three of the com- 
pounds while the fourth, diethylamine, agrees within 24 percent. The literature 
value for the critical volume of diethylamine was then used to recalculate the 
diffusion coefficient yielding a value of 0.069 cm2/sec. 


3.3 (CU) Viscosity of Vapor. 


(U) The modified Sutherland's equation shown below was developed by e 
Licht and Stechert!0 and is used to calculate the viscosity of gases. With- 

in the pressure range where viscosity is not dependent on pressure, the 

equation allows the viscosity of a gas to be calculated from its critical 

temperature, its critical pressure, and its molecular weight. 


3, 4\t/6 
a MP. o3/2 
nu = 6.30 X 10 c (6) 

where up = viscosity of vapor in cP at T°K 

M = molecular weight 

P. = critical pressure in atm 

T, = critical temperature in °K 

T. = reduced temperature = T/T, 
(U) Licht and Stechert calculated the viscosity of gases for 23 com- 


pounds and found that for 19 of the compounds, estimated values fell within 
10 percent of the literature values. The remaining compounds were water, 
helium, ammonia and ethanol. Estimated values for these compounds were 
within 21 percent of the literature values. 


(U) Table 4 lists the estimated and literature vapor viscosity values 

for the four test compounds. Using the previously estimated critical pressures 
and temperatures, the physical property data system was used to estimate the 
viscosity of vapor for these compounds. While agreement within 4 percent was 
obtained for the first three compounds, the estimated value for 2-propanol is 
20 percent lower than the literature value. The deviation was reduced to 16 
percent when the calculation was repeated using literature values for 

critical properties. These results for 2-propanol are consistent with Licht 
and Stechert's results for ethanol. ® 


(U) It should be noted that highly associated compounds may not give 
estimates as accurate as those given by the compounds found in table 4; the 
data and compounds shown in table 4 were probably used to develop the esti- 
mation techniques. 
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LIST OF COMPOUNDS (U) 


Code Number 


Code Name Name Structure 
(U) EA 1033 Bis(2-chloroethy]) sulfide; C1CH,SC,H,C1 
(HD) Distilled Mustard 
(U) EA 1034 2-Chlorovinyldichloroarsine C1CH=CHAsC1, 
(L) 
(U) EA 1036 Bis [2-(2-chloroethylmercapto) (C1CH,SC,H,) 0 
(T) ethyl] ether aR eee 
(U) EA 1053 Tris(2-chloroethy]) amine; (C1C,H,),N 
(NHS) Nitrogen Mustard -3 
(U) EA 1205S Ethyl N,N-dimethylphosphoramido- CH, 
(GA) cyanidate as 
-P-0-C_H 
cH’ |! 25 
3 CN 
(U) EA 1207 Ethyl methylphosphonofluoridate i 
GH, -F-0C 2H. 
CH 
(U) EA 1208 2-Propyl] methylphosphono- i j 3 
(GB) fluoridate CH,- aaa 
F CH, 
(U) BA 1209 2-Propyl ethylphosphono- f CHs 
(GE) fluoridate re a a 
F CH, 
i 
(U) BA 1210 Pinacolyl methylphosphono- CH,- -0-Cll-C(CH,) 3 
(GD) fluoridate 
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(U) EA 1211 2- (4-Methylpenty1)methy1- CH,-R-O-CH 
(GH) phosphonof luoridate 7 MS | CH 
CH,CH 
CH, 
(U) BA 1212 Cyclohexyl methylphosphono- afo< 
(GF) fluoridate : 
boost | 
(VU) EA 1213 2,2-Dimethylpropyl CH,- -OCH,C — CH, 
methylphosr*onof luoridate . ; eo 
3 
(U) BA 1214 2-Ethylhexy1 methy lphosphono- i cH oHs 
fluoridate CH.-F-OCH.CH 
| \ (cu) ,cH 
F 3 
i 
(U) EBA 1230 2-Phenoxyethyl methyl- CH,-P-0-C,H,-0 
phosphonof luoridate | 
F 
(U) BA 1232 Methyl methylphosphonofluoridate = - CH,-2'-0CH, 
(U) BA 1244 2-Chloroethy] N,N-dimethyl- CH 
phosphoramidocyanidate mG: -0-C.H.Cl 
“ 24 
CH 
3 CN 
(U) BA 1245 Methyl N,N-diethylphosphoramido- C,H 
cyanidate ft ccs 
f'| : e| 
CjHe CN | 
(U) BA 1246 Cyclopenty] methylphcsphono- CH,- 


fluoridate [ \ XJ 
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: (u) EA 1274 2-Pentyl methylphos - CH; - te 
| phonofluoridate No. ae CH, CH, 
) 
* (Cc) EA 1356 (Racemic) 2-methylcyclo- 2 
‘ hexyl methyl phosphono- a 
| fluoridate 
; 

(U) EA 1508 0,0-Diethyl S-(2-diethyl- CoH.0 

(VG) aminoethyl) phosphoro- > Msc gh,N(C,H,) 5 , 
thioate ete Hc a 


(U) BA 1511 3-pyridy] 3,3,S-trimethy)- ss 
(VP) cyclohexyl] methyl phosphonate 


(U) EA 1517 0-Ethyl S-(2-diethylamino- CH 
(VE) ethyl) ethylphosphonothioate i 
~SC4H,N(CjH.) > 
Cog 
oe i 
(U) EA 1576 2 -Ethoxycarbony] -1-methy] - CH, -P-0- =CH-C-0C 5H. 
vinyl 3-methylcyclohexyl 
methy! phosphonate 
0) 
Be 
a 
(U) EA 1622 O-Isopropyl S-(2-diethyl - CH, 
aminoethyl) methyl phosphono- NSP -SCQH,N(C2H.) 5 


thioate ‘4 8 
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(U) EA 1664 
(VM) 


(U) EA 1677 
(VS) 


(U) EA 1694 


(U; BA 1699 


(U) EA 1701 
(Vx) 


(U) EA 1724 
. (QL) 


(U) BA 1728 
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ethy]) methylphosphonothioate p86, H,N(CoHe) 9 
C.H 
25 
Q-Ethyl S-(2-diisopropy1- Fe) 
aminoethyl) ethylphosphono- -SC ott .s foncca) 2 
thioate CoH a : 
255 
O-Bthy] $-(2-dimethylamino- CoH | 
ethyl) ethylphosphonothioate ~SC ot,N (CH,) 5 
CHO 
O-Ethyl S-(2-dimethylamino- CH, 1 
ethyl) methylphosphonothioate NP -SC oH gN(CH3) 2 
CHO 
O-Bthyl S-(2-diisopropyl- cH, f 
aminoethyl) methyl phosphono- ~SC5H,N [owen | 2 
thioate C.H.O 
25 
0-(2-Diisopropylaminoethy]) CH 
O-ethyl methylphosphonite ' *p-00,H,N CH(CH,) >| 5 
C,H,0 
Q-Isopropyl S-(2-diisopropyl- CH, J 
aminoethyl) methylphosphono- -SC5H, N CH(CH,) 5 
thioate (Ch, ) 2cho 
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USE OF THE ANTOINE EQUATION TO FIT VAPOR PRESSURE DATA (CU) 


Vapor pressure data playa an important role in the developaent of 


critical role are listed in Table C-1l. 
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Table C-l. (U) Areas in Which Vapor Pressures Play « 
Critical Role in Chemical Defense 
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The Antoine equation! fa defined as follows: 


Log, 4? = A-B/(t+c) 


P = vwapor pressure, 
t ™ temperature in centigrade, and 


» & © = constants. 


to 230. 
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ia used extensively in this report. In xefereace 2, ft ie se 
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Cu) Table C-2 lists th. seasons why the Antoine equation hog heen 
chosen for correlating and cxtrapolating vapor pressure data. Some problem 
dealing with fitting data to the Antoine equation do eccur with the current 
computer program; ? see table.C-3. However, studies are being cenducted on 
preblems and preliminary results have been reported by Celmias.4 


Table C-2. (U) Advantages of Using the Antoine 
Equation for Extrapolations 
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Numerous comparisons with other equations have been re~ 
ported in published literature. 
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The Antoine equation: 


accurately fits the data cf hundreds of compounds over 
a wide range of temperatures. 


¢ is close to the Clevsius-Clapeyron equation. 
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Tible C-3. (U) Problems in Using the Antoine Equation 


Hl. Current computer program gives an estimate of the quality 
5 of fit, but does not give an estimate of error of specific 
vapor pregsures. 


Ze While the program fits the data very well, there are sone 
logical problems in the routine (not in the Antoine eque ion) 
that might reduce the quality of extrapolations. : 


Antoine equation extrapolations close to the critical point 
tend to be low. 
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(U3 The Antoine equation was chosen for use at CSL ater a careful ct 
of the literature and in~house testing of many equations. A number of phys 
chemists participated in the decision. A more recent study tends ¢ 

the correctness of this decision. Trump? states, “In each case the Antoine 
equating gave a cleser fit... and was therefore the preferred correlating 


‘aquation. The Chebyshev equation produced deviations about twice as large. 


ci Tne consensus® in the chemicai literature is that the Ancoine 
equation gives goad fitea at low temperaturea; but, the extrapolations up 
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reduced temperature of 0.75 with the Antoine equation result in low 
vweiues. Therefore, equations 1 and 2 (reference 2) which give even lower 
a cannot provide superior extrapolations as reported. The errors in- 
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that equation 1 (reference 2) dues not pass through any of the data points 
and that equation 2 passes through only the top two points, an insufficient 
number for the extrapolation of a curving line. Equation 2 actually goes 
through a maximum vapor pressure at the highest temperatures. 
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Figure C-1. (U) Comparison of Antoine Fquation With Other Equations. 
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